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Spectroscopically- and spatially-resolved optical line emission in the 
Superantennae (IRAS 19254-7245)* 
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ABSTRACT 

We present VIMOS integral field spectroscopic observations of the ultraluminous infrared 
galaxy (ULIRG) pair IRAS 19254-7245 (the Superantennae). We resolve Ha, [N II], [O I], 
and [S II] emission both spatially and spectroscopically and separate the emission into mul- 
tiple velocity components. We identify spectral line emission characteristic of star formation 
associated with both galaxies, broad spectral line emission from the nucleus of the southern 
progenitor, and potential outflows with shock-excited spectral features near both nuclei. We 
estimate that < 10% of the 24 fim flux density originates from star formation, implying that 
most of the 24 /am emission originates from the AGN in the southern nucleus. We also mea- 
sure a gas consumption time of ~ 1 Gyr, which is consistent with other measurements of 
ULIRGs. 

Key words: galaxies: individual: IRAS 19254-7245, galaxies: active, galaxies: ISM, galaxies: 
starburst, infrared: galaxies 



1 INTRODUCTION 

We are undertaking an integral field spectroscopic survey of 
a volume-limited sample of 18 ultraluminous infrared galaxies 
(ULIRGs; galaxies where Lir > 10 12 L(T)) with the V isible Multi- 
Object Spectrograph (VIMOS; iLeFevre et all l2003h at the Very 
Large Telescope. We have two main scientific goals. First, we want 
to study gas inflow and outflow in ULIRGs, with particular empha- 
sis on the detection of starburst-driven superwinds or AGN-driven 
jets, and determine whether these outflows may be part of a feed- 
back mechanism that inhibits star formation. Second, we want to 
examine gas excitation mechanisms within these galaxies and de- 
termine whether AGN or star formation dominate the energetics of 
ULIRGs. The results from this survey can be applied to understand- 
ing the dynamics and rest-frame optical spectra of both nearby and 
more distant ULIRGs, especially in light of recently published inte- 
gral field spectroscopic obse rvation of infrared- and subm illimetre- 
luminous z > 1 objects (e.g. Swinbank et al. 2005, 2006). 

We present here the first results from this survey: VIMOS in- 
tegral field spectroscopic observations of IRAS 19254-7245 (the 
Superantennae). The object, which has two di stinct nuclei and tw o 
tidal tails that extend over hundreds of kpc (Mirabel et_alj ll99ll) . 
has been the subject of many detailed studies. The southern nucleus 
contains an AGN, as shown by studie s at multiple wavelengths 



(iVanzi et all 120021: ICharmandaris et al.J 120021 : iBraito et al.1 120031 : 



m 

Q2D, 



Bert a et al.ll2003l :lReunanen e t alJl2007h . but the nort hern nucleus 
appears to be completely dominated by star formation dBraito et"al] 



* Based on observations collected at the European Organisation for Astro- 
nomical Research in the Southern Hemisphere, Chile [080.B-0085]. 



2003: iBerta et all [2003). Some evidence has been siven for the 
presence of several distinct dynamical components within the in- 
ner arcmin of the galaxy, including a broad line component (full- 
width half-maximum (FWHM) - 2000-2500 kins" 1 ) associ- 
ated with the southern nucleus, narrower components (FWHM < 
500 km s _1 ) as sociated with the progenitor discs, and some high- 
velocity clouds JVanzi et al.l2002l ; lReunanen et al .12 007). However, 
the previous optical studies have been mainly single slit spectra that 
have lacked spatial information, and previous near-infrared integral 
field spectroscopic observations covered only an 8 x 8 arcsec re- 
gion that does not include all the emission from the southern disc. 
The VIMOS data that we present here, which covers the central 
35 x 35 arcsec, allows us to resolve optical line emission both spa- 
tially and spectroscopically, so we can clearly discern the AGN, 
star-forming regions, and other structures within both galaxies. We 
used this galaxy to test many of the analysis methods that will be 
applied to the sample as a whole. The spectral line emission in 
this galaxy pair was detected at significantly larger radii and the 
line-emitting structures are more complex than most other sample 
galaxies, which is why we have focused on it for this first paper. 



2 OBSERVATIONS, DATA REDUCTION 

Observations were performed with the VIMOS integral field unit at 
the Very Large Telescope on 2007 October 9. The data were taken 
using the HR red grism (6300-8700 A), which has a spectral reso- 
lution of 3100, and a plate scale of 0.67 arcsec for each fibre that 
gave a field of view of 27 x 27 arcsec for each pointing. The seeing 
during the observations was ~ 0.7 arcsec. Spectra were measured 
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Figure 1. Subsections of the spectra extracted from different regions drawn as blue error bars, with the best fitting models overplotted as solid black lines. 
When two velocity components are fit to the data, component 1 is shown in red and component 2 is shown in green. For region SN, the individual spectral 
lines for the broader velocity component are shown as thin green lines. Line centers for the [O I] 6300 and 6364 A lines, the [N II] 6548 and 6583 A lines, the 
6563 A Ha line, and the 6716 and 6731 A lines are overplotted as grey lines. Parameters from these fits are given in Table [T] The spectra for region SN are 
shown over a broader wavelength range because the lines are broader than in other regions. 



Table 1. Parameters for line fits for regions in Figure [3] 



Region 


Ha 
velocity" 
(kms- 1 ) 


Ha 
line width b 
(kms" 1 ) 


[O I] 

6300 A 


Flux (erg cm 2 s 
Ha [N II] 
6563 A 6583 A 


') 

[S II] 
6716 A 


[S II] 
6731 A 


SN 


-7.9 ±0.3 
-162.9 ± 1.7 


379.8 ±0.9 
1851.3 ± 1.5 


5.4 ±0.3 
107.3 ±0.25 


67.8 ±0.3 
59 ±3 


60.0 ±0.3 
761.9 ± 1.5 


16.27 ±0.16 
53.9 ±0.6 


17.32 ±0.17 
97.9 ±0.6 


SI 


-379 ± 4 
-137 ± 14 


153 ± 7 
529 ± 17 


0.08 ±0.07 
0.48 ±0.16 


2.3 ±0.2 
3.3 ±0.2 


1.23 ±0.14 
4.5 ±0.3 


0.61 ±0.09 
1.31 ±0.17 


0.34 ±0.08 
1.10 ±0.16 


S2 


-21 ±5 
-170 ± 11 


108 ±5 
358 ±9 


0.03 ±0.05 
0.74 ±0.15 


4.5 ±0.4 
5.2 ±0.4 


1.15 ±0.15 
4.0 ±0.2 


0.53 ±0.10 
1.78 ±0.16 


0.27 ±0.08 
1.44 ±0.15 


S3 


-366.3 ±0.7 


117.3 ±0.7 


0.24 ±0.08 


11.85 ±0.14 


3.75 ±0.07 


1.6 ±0.3 


1.13 ±0.19 


S4 


256 ±6 
140 ± 20 


114 ±6 
311 ± 16 


0.05 ±0.06 
0.50 ±0.15 


2.8 ±0.4 
2.2 ±0.5 


1.2 ±0.2 
2.4 ±0.3 


0.55 ±0.13 
0.9 ±0.2 


0.33 ±0.09 
0.73 ±0.15 


NN 


34 ±3 
33 ±4 


226 ± 13 
510 ± 20 


0.5 ±0.2 
1.4 ±0.4 


7.9 ± 1.4 
11 ±3 


6.3 ± 1.3 
12 ±3 


2.0 ±0.5 
4.0 ±0.9 


1.0 ±0.3 
4.4 ±0.7 


Nl 


-4 ±6 
102 ± 10 


110 ±8 

356 ±8 


0.04 ±0.13 
0.90 ±0.15 


2.1 ±0.4 

5.2 ±0.5 


1.0 ±0.2 
4.5 ±0.3 


0.33 ±0.10 
2.01 ±0.18 


0.28 ±0.09 
1.36 ±0.15 



a These are calculated relative to a central velocity of 17950 km s 1 . Negative velocities correspond to blueshifted lines. 
b These are based on the FWHM of the lines. 



in five pointings offset from each other by 5 arcsec, which provided 
redundancy within the central region and ensures the availability of 
blank sky for background measurements. This strategy is used for 
all galaxies in the ULIRG survey. 

The VIMOS data reduction pipeline produces calibrated spec- 
tra that show the spectra measured for each fibre. The data from 
each pointing are stored in four files that contain the spectra from 



the four quadrants of the integral field unit. We used these pipeline- 
processed data for our analysis, but additional processing was 
needed to transform the spectra into a data cube. First, for each 
pointing and for each quadrant, we identified fibres that measured 
background emission by integrating the spectra and using an iter- 
ative process to remove fibres with high continuum signals, which 
would indicate the presence of emission from the target. Median 
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background spectra for each quadrant and each pointing were deter- 
mined using these background fibres, and these background spectra 
were subtracted from the data. Following this, the spectra from each 
pointing were mapped into individual spectral cubes, and then the 
cubes were median combined to produce the final spectral cube. 



3 SPECTRAL LINE FITTING AND ANALYSIS 

Before measuring spectral line emission, we corrected the data for 
redshift using a velocity of ~ 17950 km s _1 , which was estimated 
empirically to be the approximate velocity of narrow-line emission 
from the southern nucleus, and we subtracted a continuum deter- 
mined by fitting a line through the data between 6100-6200 A and 
between 6800-7000 A in the rest frame. We then fit Gaussian func- 
tions to the Ha, [O I], [N II], and [S II] lines within each spectrum 
in the data cube. We forced the fits to treat the offset between ad- 
jacent spectral features (i.e. adjacent [O I] lines, adjacent [S II], or 
lines near Ha) as constants. We generally determined whether to fit 
one or two velocity components by visually inspecting the lines for 
features such as skewed line profiles or double-peaked structures 
that would indicate that two velocity components are present, but 
we also decided to fit two velocity components to the data when 
the two component fit produced a lower reduced . When we fit 
one velocity component, we used the same line widths for adja- 
cent lines. When we fit two velocity components to the data, we 
fit both simultaneously. To reduce the uncertainties in the line fits 
where two velocity components were present, we performed fits to 
all spectral features between 6200-6800 A where the offsets among 
all spectral lines were treated as constants and where the corre- 
sponding line widths for each velocity component were treated as 
equal. We also forced all line fluxes to be positive in all fits. 

Examples of spectral line fits are shown in Figure [TJ with pa- 
rameters for the best fit lines given in Table Q] Each spectrum was 
extracted from multiple spatial pixels for the analysis later in this 
section, but similar fits were applied to the spectra for single spatial 
pixels. We also used three of these regions to show in Figure [2 ex- 
amples of the robustness of the fits to the Ha line. Region S3 is an 
example or where we observed spectral line emission from a sin- 
gle velocity component; we demonstrate with this profile that the 
spectral line is fit better by a Gaussian function than by a Lorentz 
function. Region Nl is an example of where the spectral lines are 
significantly skewed and must be fit by two Gaussian functions. Re- 
gion NN is a special case found in only one location in the object 
where the lines exhibit extra kurtosis. When we fit one Gaussian 
function, we saw that the lines had both high central peaks and 
broad wings when compared to the fit, as shown by the high posi- 
tive residuals for the Ha line fit in Figure [2] The line is better fit by 
either a single Lorentz function or by two Gaussian functions with 
similar central wavelengths. However, the emission from region Nl 
immediately to the west of NN consists of skewed lines with two 
velocity components that have line widths and line ratios similar 
to those in the two Gaussian component fit to NN. This suggest 
that the line profiles for NN should be modeled as two Gaussian 
components and not one Lorentz component. 

The results from the fits to each spatial pixel were then used 
to produce maps of the parameters, as shown in Figure[3] Addition- 
ally, Figure [4] shows how the velocities and line widths vary across 
vertical and horizontal lines that were placed so as to cover loca- 
tions of interest with two velocity components. Component 1 con- 
tains almost all of the emission from locations with only a single 
line component and most of the narrower line emission in locations 




6550 6560 6560 6570 o 6560 

Wavelength (A) 



Figure 2. The top panels show spectra for the regions around the Ha line 
extracted from the regions in Figure [3] The data are shown as blue error 
bars. The best fitting model (one Gaussian function for S3 and two Gaussian 
functions for Nl and NN) are shown in black. Individual components of the 
two Gaussian function models are shown in green. Alternate single Lorentz 
functions for regions S3 and NN are shown in orange (although the orange 
lines are mostly overlapped by the black lines); alternative single Gaussian 
functions for regions Nl and NN are shown in red. The bottom panels show 
the residuals from the best fitting models in black and the residuals from the 
alternate models in red or orange, with the individual data points shown as 
squares and the typical size of the error bars shown as a blue line on the left 
side of the panels. 

with two line components. Component 2 generally contains broader 
spectral line emission. In the southeastern HII region, however, we 
found that the broader spectral line component had velocities that 
were closer to the velocities traced by the narrow line component 
in adjacent pixels and that the velocities of the narrow component 
were sharply different from component 1 in adjacent pixels, as seen 
in Figure |4] We therefore shifted the narrower line emission for 
these locations into component 2 and the broader line emission into 
component 1 . We also found six pixels along the southern or west- 
ern edge of the detected region where either a single line compo- 
nent or a narrower line component had velocities closer to that for 
the broader line component in adjacent pixels, and so the narrower 
or single line component was shifted to component 2 in these cases 
as well. 

We can see two rotating progenitors and a couple of Ha-bright 
knots in component 1 in Figure[3] In component 2, we can see broad 
(FWHM - 1850 km s" 1 ) spectral line emission from the nucleus, 
a broad arc to the southeast of the nucleus, and a few other regions 
in both progenitors. The velocity maps and Figure HI clearly show 
that the two progenitors are rotating in the same direction. 

The individual spectra shown in Figure [T] are useful for un- 
derstanding the nature of the different regions of the galaxy. In 
the followi ng discussion , we u se the spectral line diagnostics pre- 
sented by Kew lev et aO {2006) to determine whether the spec- 
tral line emission is consistent with star formation or AGN-like 
emission, although the usefulness of these diagnostics is limited 
without additional measurements of H/3 and [O III] (5007 A) 
fluxes. We treat locations where log([Ol] (6300 A)/Ha) < -1.0, 
log([NH] (6583 A)/Ha) < and log([SlI] (6716,31 A)/Ha) < - 
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Figure 3. Images of the Ha; flux, the [N II]/Ha flux ratio, the Hen velocity, and the Hen line widths for the inner 20 x 26 arcsec of the Superantennae. The 
kinematic structures have been divided into two components, with the broader spectral lines generally placed in the second component; see the text for details. 
The objects in each frame may not necessarily be physically associated. North is up and east is left in each image, and the angular scale is shown in the Ha 
flux image for the first component. The magenta regions in the velocity images are the locations for which spectra are displayed in Figure [T] The red dotted 
lines show the locations with velocities and line widths plotted in Figure [4] 



0.2 as being dominated by photoionisation from star formation, 
whereas regions with higher line ratios are treated as having AGN- 
like emission. The southern nucleus (SN) contains both broader and 
narrower emission line components. The narrow component prob- 
ably originates from the progenitor disc remnant, and the spectral 
line ratios are consistent with star formation. The broader compo- 
nent is consistent with AGN activity. In particular, the Ha emis- 
sion is very weak compared to other spectral features, most notably 
seen in the map of the [N II]/Ha: ratio, indicating that the gas is 
ionised by a hard radiation fie l d. Thi s is consistent with results for 
Paa found bv lReunanen et al.l fa)07h . SI contains both a narrower 
component associated with the disc remnant and a broader com- 
ponent that is part of the arc extending from the southern nucleus. 
The line ratios of the broader component are consistent with AGN- 
like emission, which implies that this extension is a shock si milar 
to those seen in other ULIRGs bv lMonreal-Ibero et al.l d20Q6h . The 
apparent physical connection to the southern nucleus, the smooth 
changes in velocity from the southern nucleus to the arc, and the 
uniform [N II]/Ha ratio across the arc imply that it is gas ejected 
from the southern nucleus. However, the orientation and motion 



of the arc relative to the rotation of the southern disc remnant im- 
ply that it could be a tidal feature. The broader component in S2 
is associated with the southern disc remnant, while the narrower 
component is associated with a cloud moving at a higher velocity. 
Both components in region S2 have line ratios that are consistent 
with photoionisation. The spectrum of the region in S3 is consistent 
with star formation. The spectrum of S4 contains two components, 
and the broader component is consistent with AGN emission. In 
NN and Nl, the line ratios of the narrower components are like 
HII emission, but the broader components are AGN-like. While the 
broader components of NN and Nl may be associated with ejecta 
from the southern nucleus or with material stripped or ejected from 
the southern disc remnant, its location and velocity imply that the 
emission originates from gas ejected from the northern disc rem- 
nant. 

Based on the application of the diagnostics from lKewlev et al.l 
(2006), component 1 and the southeastern HII region in component 
2 (outlined in magenta in Figure [3]) appear to trace all of the star 
formation. The total Ha flux from star formation as traced by these 
components is (3.52 ± 0.02) x 10~ 14 erg s" 1 cm" 2 , or - 70% 
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Figure 4. Velocities and line widths for component 1 (as red diamonds) and 
component 2 (as green squares) for the locations shown by the dotted lines 
in Figure [3] 



star-forming regions in th is galaxy is not va riable and accurately 
represented by Ay — 3. 1. iBerta et al.l d2003h found that extinction 
was lower in the outer regions of the galaxy, which could reduce 
the estimated 24 /xm flux density from star formation. 

Using the extinction-corrected Ha flux for star-forming re- 
gions, a distance of 245 Mpc (calculated using a velocity of 17950 
km s _1 and Ho of 73 km s _1 Mpc -1 ), and th e conversion of Ha 
flux to star formation rate given by iKennicuttl dl998h , we estimate 
the star formation rate to be 25 M yr _1 . The total molecular gas 
mass has been measured to be 1.9 — 3.0 x 10 10 M dMirabel et al.1 
1990: IVanzi et al.l l2002). which would imply a gas consumption 
time of ~ 1 Gyr. While this gas consumption time does not account 
for gas recycling, it is still indicative of the efficiency of star for- 
mation. This is slightly high but still consistent with the "few times 
10 8 yr" gas consumption times estimated fo r ULIRGs using sta r 
formation rates derived from far-infrared data jTacconi et al .12 006), 
but, as expected, it is lower than the ~ 3 Gyr gas consumpti on times 
measured for normal spiral galaxies dKennicutt et al .11 1994b . 



of the total Ha flux from the galaxy. To understand the relative 
contribution of star formation to the 24 /xm flux density, we can use 
extinction measurements from the literature and the Ha flux in 

fH*e°- 812Av = f Ha + 0.022-^- f 24fim (1) 

24/xra 

(adapted from IZhu et al.l d2008h ) to estimate the 24 /xm flux den- 
sity from star formation. In this equation, fua is the measured 
Ha flux, fHc<e ' 812Av represents the extinction-corrected Ha flux, 
and /24/im is the 24 /im flux density. Although similar equations 
have been publi shed for individual HII regions within galaxies, the 
IZhu et al.1 ( 2008) version is the only one currently published that 
has been calibrated for use on global flux measurements, and given 
that extreme starbursts and AGN with L/i?uptolO 12 L were in- 
cluded in their analysis, their relati on should be applicable t o this 
galaxy. The extinction function of Savage & Mathisl dl979h with 
R = 3.1 was used to derive the extinction correction term on the 
left side of the equation. The Ay measured in the southern nu- 
cleus bv lVanzi et all d2002h and lBerta etaD (120031) using single-slit 
spectra is ~ 3.1. Assuming that this extinction applies to all star- 
forming regions and that the extinction is intrinsic to the source 
itself, we estimate the 24 /xm flux density from star-forming re- 
gions to be 0.15 Jy. In contrast, the total 25 /xm flu x density 
measured by IRAS is 1.24 ± 0.06 Jy dMoshir et all 19891) . We ten- 
tatively conclude that >90% of the 24 /im flux originates from 
something other than star formation; the AGN in the southern nu- 
cleus is the most likely source of this emission. 

Our results are cons istent with the spe ctral energy distribu- 
tion (SED) modelling by IBerta et al.l (120031) . which predicts that 
the AGN is the dominant source of 25 /xm emission in the Super- 
antennae. While iFarrah et al.l d2003h found using SED template fit- 
ting that AGN may be the primary source of 25 /im emission in 
ULIRGs in general, we disagree with their conclusions that the Su- 
perantennae SED can be explained entirely by star formation. The 
starburst template that they used does not accurately describe the 
Superantennae data between 10-100 /xm, so it may be an inaccurate 
description of the 24 /xm flux density. While the 24 /xm flux den- 
sity may be dominated by an AGN, star formation may represent 
approximately half of 1-1000 /xm flux, and that starburst emission 
may dominate the far-infrared em ission from the Superantennae, 
as predicted by Berta et al Our results rely on the assump- 

tion that the kinematic components that we have identified are the 
only locations with star formation and that the extinction across the 
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